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Using Drell-Yan pair production as an example, we explore the relation between two well-known 
mechanisms for single transverse-spin asymmetries in hard processes: twist-three quark-gluon cor- 
relations when the pair's transverse momentum is large, q± ^ Aqcd, and time-reversal-odd and 
transverse-momentum-dependent parton distributions when qx is much less than the pair's mass. 
We find that although the two mechanisms have their own domain of validity, they describe the 
same physics in the kinematic region where they overlap. This unifies the two mechanisms and 
imposes an important constraint on phenomenological studies of single spin asymmetries. 



1. Single-transverse spin asymmetries (SSAs) have a 
long history, starting from the 1970s and 1980s when 
large SSAs were first observed in hadronic reactions at 
forward rapidities. The size of the asymmetries came 
as a surprise and has posed a challenge for researchers 
in Quantum Chromodynamics (QCD). SSAs have again 
attracted much interest in recent years, both experimen- 
tally and theoretically Two mechanisms have been 
proposed in QCD for explaining the large SSAs: effects 
due to time-reversal-odd and transverse-momentum de- 
pendent parton distributions (Sivers functions) 0, or 
to spin-dependent twist-three quark-gluon correlations 
(Efremov-Teryaev-Qiu-Sterman (ETQS) mechanism) 0, 
IJ]. Much progress has been made in understand- 
ing theoretical issues involved in transverse-momentum- 
dependent (TMD) parton distributions. Their gauge 
properties have been clarified H, IE IE S and 
the relevant QCD factorization formulas have been es- 
tablished IE lEi Qi ^3 • Phenomenological studies have 
been made to explain the available data, using both 
mechanisms 0, 01 . An interesting relation between the 
two types of distribution functions was discussed in |lOl| . 
However, a clear connection between the two seemingly- 
different physical mechanisms, particularly at the level of 
physical observables, has so far not been established (see 

In this letter we explore the relation between the two 
mechanisms in Drell-Yan pair production at transverse 
momentum q^. If the pair's mass is Q, the quark-gluon 
correlation mechanism works when q_\_. Q ^ Aqcd , the 
strong interaction scale. On the other hand, the TMD 
QCD factorization formalism applies when q±^ <^ Q. 
Therefore, one expects there to be a common kinematic 
region, Aqcd ^ 9j_ ^ Q, where both mechanisms should 
apply and describe the same physics. We perform calcu- 
lations in this intermediate region in both formalisms and 
find that they indeed yield the same result. This conclu- 
sion imposes a rigorous constraint on phenomenological 
studies of the single-spin asymmetry data. 



Consider scattering of a transversely polarized proton 
of spin Sa_ and momentum P on an unpolarized hadron 
(another proton, for example) of momentum P' , produc- 
ing a virtual photon that subsequently decays into a pair 
of leptons. The total center-of-mass energy is denoted 
by s = (P + P')2. At the lowest order in QCD per- 
turbation theory, there are two ways to produce such 
a pair: either a quark-antiquark pair annihilates into a 
virtual photon after emitting a gluon of transverse mo- 
mentum q_\_, or a quark scatters with a gluon producing 
a quark and a virtual photon. If the lepton pair has a 
large positive rapidity, y (in the forward direction of the 
polarized proton) , the dominant contribution comes from 
scattering involving a quark from the polarized proton. 
Our interest is to calculate the spin- dependent differen- 
tial cross section dA(j{S±)/dQ^dyd^q± at large y with 

2. When q±, Q ^ Aqcd, the spin-dependent cross 
section can be calculated in terms of a twist-three quark- 
gluon correlation. The physics of this correlation can be 
seen as follows: when a transversely-polarized proton is 
traveling at nearly the speed of light, its internal color 
electric and magnetic fields have preferred orientations 
in the transverse plane. By parity invariance, the color 
electric field must be orthogonal to the spin of the proton. 
If averaged over the proton wave function, the field van- 
ishes because the proton is color-neutral (also because 
of time-reversal symmetry). However, if one multiplies 
the color electric field with the quark color current, the 
average may be non-zero. This average defines a quark- 
gluon correlation function that characterizes a property 
of a polarized proton. 

High-energy scattering probes the so-called light-cone 
correlations. For these, quark and/or gluon fields are sep- 
arated along the light-cone direction (if denotes a 
space-time coordinate, the light-cone variables are de- 
fined as C± = ie±e)/V2,U = A typical 
diagram for scattering of an antiquark in the polarized 
gluon background is shown in Fig. 1. The blob in the 



2 























P.S_ 





FIG. 1: A generic Feynman diagram contributing to the single 
transverse-spin asymmetry for the Drell- Yan process through 
the quark- antiquark scattering channel, containing both soft- 
and hard-pole contributions. 



lower part of the diagram represents the spin-dependent 
quark-gluon correlation described above. One finds the 
following expression for the correlation 0: 



Tf{xi,X2) 



c?C dr] 
47r 



(1) 



where the sums over color and spin indices are implicit, 
\PS) denotes the proton state, ip the quark field, and 
the gluon field tensor. In Eq. xi ~ k'^^/P^ 
and X2 — fc^/-P^ are the fractions of the polarized pro- 
ton's light-cone momentum carried by the quark in Fig. 1, 



while 



= fc+/P+ - X2 



Xl 



is the fractional mo- 



mentum carried by the gluon; C is the light-cone gauge 

link, £(^2,Ci) — 6xp ^— that makes 

the correlation operator gauge- invariant; and e°f^ is the 
2-dimensional Levi-Civita tensor with — 1. 

When the antiquark of the unpolarized proton scat- 
ters off a quark of the polarized proton in the presence 
of a polarized color electromagnetic field (the gluon) in 
Fig. 1, the lepton-pair yield has a dependence on the 
transverse orientation of the spin. The strong interac- 
tion phase necessary for SSAs arises from the interfer- 
ence between an imaginary part of the partonic scatter- 
ing amplitude with the extra gluon and a real part of 
the scattering amplitude without a gluon in Fig. 1. The 
imaginary part is from the pole of the parton propagator 
associated with the integration of the gluon momentum 
Xg. Depending on which propagator's pole contributes, 
Aa{S±) gets contributions from Xg — (soft-pole) 
and Xg ^0 (hard-pole) 

3. By considering all partonic diagrams in Fig. 1 
and their complex conjugates, one finds for the spin- 
dependent differential cross section: 

dQ"'dya^q± J x x 

xC-[H^ + H^]ix,s,i,u) , (2) 

where the sum over all quark flavors, weighted with 
their electric charge squared, is implicit. We have de- 
fined the factor C = age"^ S±aq±i3C(s/2TT'^ with (Tq ~ 



4:TTa'^^^/3NcsQ'^ , Nc = 3. x' is the momentum fraction 
carried by the parton from the unpolarized proton, and 
X is the total parton momentum fraction from the po- 
larized proton while Xg is fixed using the pole (on-shell) 
condition. The partonic Mandelstam variables are de- 
fined as s = {xP + x'P')^J^ [xP~qf,u^ {x'P'-qf. 

and are the soft- and hard-pole contributions, 
respectively. 

The soft-pole contributions can be calculated in the 
same way as for direct-photon and inclusive-hadron pro- 
duction 4]. A total of eight diagrams sum to 
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-Tf{x,x)^ 
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(3) 



where the hard coefhcients D^^ and N^^ are 



r^<! u t 

D - = - + - 



2QH 
ut 



gg 



1 



[Q\u^ - P) + 2Q2j(g2 _ 2i) 



with l/2Nc the color factor. In the real photon limit, 
= 0, we obtain the direct-photon cross section. 
Kinematics allows only a i-channel propagator to con- 
tribute to a hard pole, and there are a total of 12 such 
diagrams |l6|. The sum of their contributions is 



H^=TF{x,x-Xg 



1 



Cf- 



(4) 



where 



qq 



[Q^ - if + 1 



As mentioned above, for the hard-pole contributions 
Xg — —xt/{Q'^ — i) differs from zero. We find that 
the hard-pole contribution has no "derivative-term" , oc 
OTf^x, x) I dx^ due to a cancellation among the diagrams. 

4. When q±^ <^ Q, single spin asymmetries can 
be generated from a spin-dependent TMD quark dis- 
tribution introduced by Sivers 0- The physics of the 
Sivers function can be understood as follows. Consider 
a transversely-polarized proton traveling with large mo- 
mentum P. The distribution of quarks with longitudinal 
and transverse momentum xP and k± can have a depen- 
dence on the orientation of k± relative to the polariza- 
tion vector S±. The TMD quark distributions are best 
defined through the following matrix: 
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(2^)2 



X (P .Z'f(r,0,6^)7+'Z'„"(0) P 



(5) 



Here the vector p = (p^,0 ,0^) is along the 
momentum direction of the proton, tf'u(C) = 




(a) (6) (c) 

FIG. 2: Decomposition of the generic diagram in Fig. 1 into 
different regions in the transverse-momentum-dependent fac- 
torization approach. The large transverse momentum of the 
Drell-Yan pair may come from (a) the Sivers function, (h) the 
antiquark TMD distribution, and (c) the soft factor. 
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FIG. 3: A typical Feynman diagram contributing to the Sivers 
functions at a large transverse momentum, k±_ ^ Aqcd, aris- 
ing from the twist-three quark-gluon correlation. The bars in- 
dicate the soft and hard poles. 



oo; and the gauge hnk is £„(— cx);^) = 

exp (^—ig /q °° dXv ■ A{\v + . In a system where 

P+ ^ , the vector v with ^ ,v±_ = is 
given a non-zero component in order to regulate the 
light-cone singularities. An expansion of the matrix in 
Dirac indices yields 

X = i k^) i, + qrix, k^)e^,^prp''k''S'' + ...] 

_ (6) 

where q{x,k±) is the spin-independent distribution and 
qrix, k±) is the Sivers function. The dependence on v is 
of the form = {2v ■ P)'^ /v"^ and is implicit. 

When q± ^ Q, the Drell-Yan cross section in lead- 
ing order q±/Q can be factorized in terms of TMD quark 
distributions. A generic Feynman diagram can be decom- 
posed into various contributions, each of which is repre- 
sented by one of the reduced diagrams shown in Fig. 2 
0. This factorization property can be seen from the di- 
agram in Fig. 1 by considering the real gluon emitted 
either along the direction of P or P'. The single-spin- 
dependent part of the Drell-Yan cross section has the 
following factorization 

d'A<j{S^) e-^S^^q^f, f 2^ 2r rr f.^ 

dQ^dyd^qr ""' A/P J ^ '^^^ ^^"^ 
,,h^5^^\k,^ + k^^ + \^-q^) 



X 9t (2^1,^1 J., Ci) 9(2:2,^2^,(2) S{X^) 



-1 



where zi = Q/y^e'^ and Z2 — Qj 'Jse~'^ are the momen- 
tum fractions of the colliding hadrons associated with 
the observed Drell-Yan pair. 7? is a hard factor and is 
entirely perturbative. The soft-factor S' is a vacuum ma- 
trix element of Wilson lines and captures the effect of 
soft gluon radiation. Since the soft-gluon contribution in 
the TMD distributions has not been subtracted, the soft- 
factor comes with inverse power. The soft-gluon rapidity 
cut-off p is defined s& p — \J (2v\ ■ v-i)^ lv\v\, where v\ 
and vi are the directions of the gauge links in the Sivers 
and the unpolarized antiquark distributions. In a special 
coordinate frame, z\C^ = z\C^ = pQ'^ 



5. From what we have discussed so far, the two mecha- 
nisms for the single-spin asymmetries appear to be quite 
different. In the case of the twist-three correlation, the 
nonperturbative function itself is time-reversal even. On 
the other hand, in the TMD factorization, the Sivers 
function is time-reversal odd and inherently contains the 
final-state interaction through a Wilson line. There is, 
however, a common kinematic region, Aqcd ^ <C Q, 
where both mechanisms should work. In this region, q±^ 
is large enough so that the asymmetry is a twist-three 
effect. But at the same time, q±^ ^ Q, so the TMD 
factorization formalism also applies in this region. 

In order to see the consistency of the two approaches, 
we need to calculate the explicit gj^-dependence in the 
TMD factorization. To do that, we let one of the trans- 
verse momenta fc^^ and be of the order of q±^ and the 
others small. When A^ is large, for example, we neglect 
ki±^ in the delta function, and the integrations over these 
momenta yield either the ordinary antiquark distribution 
or a moment of the Sivers function. The latter is related 
to the twist-three correlation lOl: 



d^kj_q{x,k±_) = q{x) , / d^k^—j^qriki^.x) ^Tf{x,x) 

iVlp 



When one of the ki± is taken to be of order and 
A^ is neglected in the delta function, one also needs 
/ d'^\A^S{\±) = 1. 

What remains is to calculate the large-fc_L behavior of 
the TMD distributions and the soft factor. The latter 
has been determined in Ref. flU and reads 

S{\^) = <5(2)(A^) + ^^Cf (lnp2 - 2) , (8) 

which is infra-red finite. To the same order, the unpolar- 
ized antiquark TMD distribution is given by 



q{z2,k^) = q{z2)5^^\k^) + ^^Cf 



(9) 



1 + 



2) + 



I) In 



kl 



where ^2 = ^2/2^' and an additional contribution involv- 
ing the gluon distribution has been neglected (IQ] . 
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As discussed above, the Sivers function at large k± 
can be calculated in terms of the twist-three quark-gluon 
correlation. There are also soft and hard pole contribu- 
tions in this case. We show one example of the relevant 
Feynman diagrams in Fig. 3. There are a total of eight 
diagrams for the soft-pole contribution, and twelve for 
the hard-pole one. Adding all contributions, wc find: 



2Mp 



47r2 (fc2^)2 



X Si^i - 1) In 



J ^{A + CfTp{x,x) 



(10) 



where the l/(/c^)2 behavior follows from a dimensional 
analysis, and where 



A-- 



1 



d 



2N, 

+Tf{x,x) 



g^TFix,x)il + + Tf{x, X - ^g)-(J-ZJ^ 
(1-6)^(26 + 1)- 



+CfTf{x, X ~ Xg) 
with X 



1 + 6 



(i-a)H 



(11) 



3 - (1 - S,i)x and 6 = zi/x. 
Plugging the above results into the TMD factorization 
formula ((7|), we obtain the spin-dependent cross section 
for the q + q channel at large q± : 



dQ'^dy(Pq± 



C 



— ^q{x') 



(12) 



X {5(6- 1)^ + 5(6-1)5} , 



where A and C have been given above and 



B — CfTf{x, x) 



1 



CI 



(1-6)+ 



+ 25(6 - 1) In 



where 6 = z^j x' . To compare this to the twist-three cor- 
relation calculation, we expand Eq. (2) for small qs^/Q. 
At leading order in q±/Q, the partonic Mandelstam vari- 
ables become s = ql/{l - 6)(1 - 6), i = -9i/(l - 6), 
and u = — gj/ (1 — 6)- Using the above, we find that the 
spin-dependent cross section to leading order in q±/Q is 
exactly the same as the one shown in (|12|l . The same 
conclusion applies for the contribution to the single-spin 
asymmetry generated by gluon-quark scattering. 

To summarize, we have studied the single transverse- 
spin asymmetry in Drell-Yan pair production at both 
large and small transverse momenta q± of the lepton pair. 
At large q±, one has a result in terms of a twist-three 
quark-gluon correlation function. At small q±, a factor- 
ization formula in terms of TMD parton distributions is 
valid. We have demonstrated that in the intermediate 
region, both approaches give the same answer. This es- 
tablishes in an explicit calculation for a physical process 
the connection between these two mechanisms for gener- 
ating single-spin asymmetries. It is this connection that 



unifies the physical pictures for the underlying dynam- 
ics of single transverse-spin asymmetries and imposes an 
important constraint on phenomenological studies of the 
existing and future data. Extension to the gluon initi- 
ated subprocess and other physical processes, such as the 
semi-inclusive deep inelastic scattering, will be presented 
elsewhere, along with the details of the calculations de- 
scribed above. 
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